Structural relaxation due to electronic correlations in the paramagnetic insulator 
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A computational scheme for the investigation of complex materials with strongly interacting 
electrons is formulated which is able to treat atomic displacements, and hence structural relaxation, 
caused by electronic correlations. It combines ab initio band structure and dynamical mean-field 
theory and is implemented in terms of plane-wave pseudopotentials. The equilibrium Jahn- Teller 
distortion and antiferro-orbital order found for paramagnetic KCuFs agree well with experiment. 
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In materials with correlated electrons, the interaction 
between spin, charge, orbital, and lattice degrees of free- 
dom leads to a wealth of ordering phenomena and com- 
plex phases [l|. The diverse properties of such systems 
and their great sensitivity with respect to changes of ex- 
ternal parameters such as temperature, pressure, mag- 
netic field, or doping also make them highly attractive 
for technological applications [l|. In particular, orbital 
degeneracy is an important and often inevitable cause for 
this complexity [2]. A fascinating example is the coop- 
erative Jahn- Teller (JT) effect — the spontaneous lifting 
of the degeneracy of an orbital state — leading to an oc- 
cupation of particular orbitals ("orbital ordering") and, 
simultaneously, to a structural relaxation with symmetry 
reduction. 

The electronic structure of materials can often be de- 
scribed quite accurately by density functional theory in 
the local density approximation (LDA) [Sl or the gen- 
eralized gradient approximation (GGA) [J, |5|]. How- 
ever, these methods usually fail to predict the correct 
electronic and structural properties of materials where 
electronic correlations play a role. Extensions of LDA, 
e.g., LDA-I-U and SIC-LDA 01 can improve the re- 
sults, e.g., the band gap value and local moment, but 
only for systems with long-range order. Hence the com- 
putation of electronic, magnetic, and structural prop- 
erties of strongly correlated paramagnetic materials re- 
mains a great challenge. Here the recently developed 
combination of band structure approaches and dynamical 
mean-field theory [8|, the so-called LDA-I-DMFT com- 
putational scheme |9[, has become a powerful new tool 
for the investigation of strongly correlated compounds in 
both their paramagnetic and magnetically ordered states. 
This technique has recently provided important insights 
into the properties of correlated electron materials [lO|, 
especially in the vicinity of a Mott metal-insulator tran- 
sition as encountered in transition metal oxides [l| . 



Applications of LDA+DMFT so far mainly employed 
linearized and higher-order muffin-tin orbital [L(N)MTO] 
methods [ill ] and concentrated on the study of corre- 
lation effects within the electronic system for a given 
ionic lattice. On the other hand, the interaction of the 
electrons with the ions also affects the lattice structure. 
LDA-I-DMFT investigations of particularly drastic exam- 
ples, the volume collapse in paramagnetic Ce V^ ll3| and 
Pu [iJI and the magnetic moment collapse in MnO [15j . 
incorporated the lattice by calculating the total energy 
of the correlated material as a function of the atomic 
volume. However, for investigations going beyond equi- 
librium volume calculations, e.g., of the cooperative JT 
effect and other subtle structural relaxation effects, the 
L(N)MTO method is not suitable since it cannot deter- 
mine atomic displacements reliably. This is partly due 
to the fact that the atomic-sphere approximation used in 
the L(N)MTO scheme, with a spherical potential inside 
the atomic sphere, completely neglects multipole contri- 
butions to the electrostatic energy originating from the 
distorted charge density distribution around the atoms. 
By contrast, the plane- wave pseudopotential approach 
employed here does not neglect such contributions and 
can thus fully describe the effect of the distortion on the 
electrostatic energy. 

In this Letter, we present a computational scheme 
which allows us to calculate lattice relaxation effects 
caused by electronic correlations. To this end, the 
GGA+DMFT — a merger of the GGA and DMFT — 
is formulated within a plane-wave pseudopotential ap- 
proach [16|, [iTI, |l8[. Thereby the limitations of the 



L(N)MTO scheme in the direct calculation of total en- 
ergies are overcome. In particular, we apply this new 
method to determine the orbital order and the coopera- 
tive JT distortion in the paramagnetic phase of the pro- 
totypical JT system KCuFa. 

KCuFs has long been known to be a prototypical mate- 



rial with a cooperative JT distortion [3| where the elec- 
tronic degrees of freedom are the driving force behind 
the orbital order [3, [g, [lM|- Indeed, the relatively high 
(tetragonal) symmetry makes KCUF3 one of the simplest 
systems to study. In particular, only a single internal 
structure parameter, the shift of the in-plane fluorine 
atom from the Cu-Cu bond center, is needed to describe 
the lattice distortion. Moreover, there is only a single 
hole in the d-shell so that complications due to multiplet 
effects do not arise. KCuFs is an insulating pseudocubic 
perovskite whose structure is related to that of high-Tc 
superconductors and colossal magnetoresistence mangan- 
ites. The copper ions have octahedral fluorine surround- 
ing and arc nominally in a Cu^+ (Sd^) electronic config- 
uration, with completely filled t2g orbitals and a single 
hole in the eg states. The cubic degeneracy of the Cu eg 
states is lifted due to a cooperative JT distortion leading 
to an elongation of the CuFg octahedra along the a and b 
axis, and an antiferro-distortive pattern in the ab plane 
[20| . This is associated with an alternating occupation of 
dx^-z^ and dy2_^2 hole orbitals along the a and b axes, re- 
sulting in a tetragonal compression {c/a < 1) of the unit 
cell. Purely electronic effects as in the Kugel-Khomskii 
theory [4I and the electron-lattice [21| interaction have 
been discussed as a possible mechanism behind the or- 
bital ordering in KCuFa. The antiferro (a- type) and fer- 
rolike (d-type) stacking of the ab planes along the c axis 
give rise to two different structural polytypes, which have 
been identified experimentally at room temperature [22| . 

Below the Ncel temperature (Tat ^38 K for a- type 
and ~22 K for d-type ordering), which is much lower 
than the critical temperature for orbital ordering, KCuFs 
shows A-type antiferromagnetic order [23|. The antifer- 
romagnetic structure is consistent with the Goodenough- 
Kanamori-Anderson rules for a superexchange interac- 
tion with dx2-z^/dy2_22 antiferro-orbital ordering. This 
is also found within LDA-I-U which finds the correct 
orbitally ordered, antiferromagnetic insulating ground 
state [y, |2J], while LDA predicts metallic behavior. 
Moreover, LDA+U calculations for a model structure 
of KCuFa in which cooperative JT distortions are com- 
pletely neglected reproduce the correct orbital order, sug- 
gesting an electronic origin of the ordering Jy, [l9| in 
agreement with the Kugcl-Khomskii theory [3. Alto- 
gether, LDA+U is able to determine the JT distortion in 
KCuFs rather well [6|, |2J| , but simultaneously predicts an 
additional long-range magnetic order. Therefore LDA+U 
cannot explain the properties at temperatures above Tjy 
and, in particular, at room temperature, where KCuFa is 
a correlated paramagnetic insulator with a robust JT dis- 
tortion which persists up to the melting temperature. To 
determine the correct orbital order and cooperative JT 
distortion for a correlated paramagnet, i.e., to perform a 
structural optimization, we here employ GGA+DMFT. 

We first calculate the GGA band structure of KCuFs 
at room temperature (space group 14/ mem) [20|, em- 



ploying the plane-wave pseudopotential approach [l7|, 
|25| . Calculations are performed for values of the in-plane 
JT distortion 6jt UM ranging from 0.2% to 7% while 
keeping the lattice parameters a and c and the space 
group symmetry fixed. In the paramagnetic phase, and 
for all values of Sjt considered here, the GGA yields a 
metallie rather than the experimentally observed insulat- 
ing behavior, with an appreciable orbital polarization due 
to the crystal field splitting. Overall, the GGA results 
qualitatively agree with previous band-structure calcula- 
tions [6|,|24|. Obviously, a JT distortion by itself, without 
the inclusion of electronic correlations in the paramag- 
netic phase, cannot explain the experimentally observed 
orbitally ordered insulating state of KCUF3. 

To include the electronic correlations, we construct an 
effective low-energy Hamiltonian Hqqa for the partially 
filled Cu Cg orbitals for each value of the distortion 5jt 
considered here. This is achieved by employing the pseu- 
dopotential plane-wave GGA results and making a pro- 
jection onto atomic-centered symmetry-constrained Cu 
eg Wannier orbitals [l6|. Taking the local Coulomb re- 
pulsion U and Hund's rule exchange J into account, one 
obtains the following low-energy Hamiltonian for the two 
(to== 1,2) 



Cu eg bands: 



H = Hgga + uJ2 
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Here the second and third terms on the right-hand side 
describe the local Coulomb interaction between Cu Cg 
electrons in the same and in different orbitals, respec- 
tively, with V = U — 2 J, and Hjjc is a double counting 
correction which accounts for the electronic interactions 
already described by the GGA (see below). To compute 
the electronic correlation induced structural relaxation of 
KCuFa, we calculate the total energy as |13l.ll8j 



E = Egga[p] + {Hgga) -^e. 
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where Egga[p] is the total energy obtained by GGA. The 
third term on the right-hand side of Eq. (2) is the sum of 
the GGA Cu Cg valence-state eigenvalues and is given by 
the thermal average of the GGA Hamiltonian with the 
GGA Green function G^'^^{iiUn): 
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{Hgga) is evaluated similarly but with the full Green 
function including the self-energy. The interaction energy 
{Hjj) is computed from the double occupancy matrix. 
The double-counting correction Eoc = \UNeg{Nf^^ — 
I) — jJNeg {Ne^ — 2) corresponds to the average Coulomb 
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FIG. 1: (colour online) Orbitally resolved Cu Cg spec- 
tral densities of paramagnetic KCuFa as obtained by 
GGA+DMFT(QMC) for different values of the JT distortion. 
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FIG. 2: (colour online) Comparison of the total en- 
ergies of paramagnetic KCuFa computed by GGA and 
GGA+DMFT(QMC) as a function of the JT distortion. Er- 
ror bars indicate the statistical error of the DMFT(QMC) 
calculations. 
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repulsion between the Ni 
orbitals. 

The many-body Hamiltonian (1) 
DMFT for t/ = 7 eV and J = 0.9 eV 
Monte Carlo (QMC) calculations [I 
shows the spectral density of paramagnetic KCuFs, ob- 
tained from the QMC data by the maximum entropy 
method, for three values of the JT distortion Sjt- Most 
importantly, a paramagnetic insulating state with a 
strong orbital polarization is obtained for all 5jt- The 
energy gap is in the range 1.5-3.5 eV, and increases with 
increasing 6jt- The sharp feature in the spectral density 
at about —3 eV corresponds to the fully occupied 3z'^~r'^ 
orbital [30[ , whereas the lower and upper Hubbard bands 
are predominantly of x^ — y^ character and are located 
at —5.5 eV and 1.8 eV, respectively. 

The total energies as a function of the JT distortion 
obtained by the GGA and GGA-hDMFT, respectively, 
are compared in Fig. [21 We note that the GGA not only 
predicts a metallic solution, but its total energy is seen to 
be almost constant for < 5jt ^ 4%. Both features are 
in contradiction to experiment since the extremely shal- 



low minimum at 6jt — 2.5% would imply that KCuFa 
has no JT distortion for T > 100 K. By contrast, the in- 
clusion of the electronic correlations among the partially 
filled Cu eg states in the GGA+DMFT approach leads to 
a very substantial lowering of the total energy by ~ 175 
meV per formula unit. This implies that the strong JT 
distortion persists up to the melting temperature (> 1000 
K), in agreement with experiment. The minimum of 
the GGA-I-DMFT total energy is located at the value 
Sjt = 4.2%, which is also in excellent agreement with 
the experimental value of 4.4% [20|. This clearly shows 
that the JT distortion in paramagnetic KCuFs is caused 
by electronic correlations. 

An analysis of the occupation matrices for the eg Cu 
Wannier states obtained by the GGA-I-DMFT calcula- 
tions confirms a substantial orbital polarization in the 
calculated paramagnetic phase of KCuFs. As shown in 
Fig. m the orbital order parameter (defined as the differ- 
ence between 3z^—r^ and x^ —y^ Cu eg Wannier occupan- 
cies [30|) saturates at about 98% for bjT ^ 4%. Thus, the 
GGA-I-DMFT result shows a predominant occupation of 
the Cu Zz^ — T^ orbitals. We note that even without a 
JT distortion the orbital order parameter would remain 
quite large (~40%). Moreover, while the GGA result 
for bjT = yields a symmetric orbital polarization with 
respect to C4 rotations around the c axis, spontaneous 
antiferro-orbital order is found in GGA-I-DMFT. This 
difference is illustrated in Fig. [31 where insets (a) and 
(c) depict the hole orbital order obtained by the GGA 
and GGA-hDMFT for Sjt = 0.2%, respectively. The 
GGA charge density is more or less the same along the a 
and b axis [inset (a)]; i.e., the Cu d^^-^i and dy2_z2 hole 
orbitals are almost equally occupied and hence are not 
ordered. By contrast, the GGA+DMFT results clearly 
show an alternating occupation [inset (c)], corresponding 
to the occupation of a x'^ — y"^ hole orbital in the local 
coordinate system, which implies antiferro-orbital order. 
For the experimentally observed value of the JT distor- 
tions of Sjt = 4.4%, both GGA and GGA+DMFT find 
antiferro-orbital order [insets (b),(d)]. However, we note 
again that, in contrast to the GGA+DMFT, the GGA 
yields a metallic solution without any JT distortion for 
T ^ 100 K, in contradiction to experiment. 

In conclusion, by formulating GGA+DMFT — the 
combination of the ah initio band structure calculation 
technique GGA with the dynamical mean-field theory — 
in terms of plane- wave pseudopotentials [l6|, we con- 
structed a robust computational scheme for the inves- 
tigation of complex materials with strong electronic in- 
teractions. Most importantly, this framework is able to 
determine the correlation induced structural relaxation 
of a solid. Results obtained for paramagnetic KCuFa, 
namely an equilibrium Jahn- Teller distortion of 4.2% and 
antiferro-orbital ordering, agree well with experiment. 
The electronic correlations were also found to be respon- 
sible for a considerable enhancement of the orbital po- 
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FIG. 3: (colour online) Dependence of the orbital order pa- 
rameter in paramagnetic KCuFs on the JT distortion as ob- 
tained by GGA and GGA-fDMFT(QMC), respectively. Er- 
ror bars indicate the statistical error of the DMFT(QMC) 
calculations. Insets (a)/(b) refer to GGA and (c)/(d) to 
GGA+DMFT results and show the hole orbital ordering for 
5jT =0.2%/4.4% (see text). 



larization. The GGA+DMFT scheme presented in this 
paper opens the way for fully microscopic investigations 
of the structural properties of strongly correlated elec- 
tron materials such as lattice instabilities observed at 
correlation-induced metal-insulator transitions. 
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